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Abstract 

In an effort to investigate some of the low energy properties of QCD, in particular 
those related to chiral symmetry breaking, as well as to obtain insights on the behavior 
of an interacting theory of fermions on the lattice, the two flavor Nambu-Jona-Lasinio 
model with SU{2) x SU (2) chiral symmetry is studied on the four-dimensional hyper- 
cubic lattice using large N techniques and numerical simulations. Naive and Wilson 
fermions are considered and transparent results are obtained regarding the following: 
the scalar and pseudoscalar spectrum, the approach to the continuum and chiral lim- 
its, the size of the corrections, and the effects of the zero momentum fermionic 
modes on finite lattices. Also, some interesting observations are made by viewing the 
model as an embedding theory of the Higgs sector. 
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When the high frequency modes of the gauge and fermionic fields of QCD are integrated 
down to the energy scale E corresponding the correlation length of the gauge field {E ^ 
glueball mass ~ 1550 MeV the resulting effective theory will essentially be a theory 
of fermions with contact interactions and cutoff A < 1550 MeV. The resulting effective 
Lagrangian will maintain the original chiral symmetry but will be more complicated. At 
energies much smaller than A it should be enough to keep in the Lagrangian the least 
irrelevant operator, namely the four-Fermi dimension six operator. This is one way to 
motivate the Nambu-Jona-Lasinio (NJL) model. 

Unfortunately, by only keeping the four-Fermi operator, valuable information was lost 
and the NJL model does not confine the quarks. Furthermore, if, for example, we want 
to study the a particle, which on phenomenological grounds is believed to have mass ~ 
750 MeV, then the separation of scales is probably not large enough to justify the neglect of 
operators with dimension higher than six. Nevertheless, the NJL model possesses the same 
chiral symmetry as QCD and it can realize this symmetry in the Goldstone mode. It is this 
feature that is most crucial in the understanding of the lightest hadrons. Furthermore, our 
interest in the model is not so much aimed at its quantitative predictability but rather on 
the qualitative insights that can provide as a low energy theory of QCD, as an embedding 
theory of the Scalar Sector of the Minimal Standard Model and as a four-dimensional 
interacting theory of fermions on the lattice. 

The NJL model has been studied extensively for various cases with continuum type 
regularizations. For a comprehensive review the reader is referred to [0] and references 
therein. The model has also been studied on the lattice in connection with the possible 
equivalence of the top quark condensate with the Higgs field [Q. In that work, however, the 
separation of scales is very large (the cutoff is of the order 10^^ GeV), and it is therefore 
quite a different problem than the one considered here. 

In this paper we consider the two flavor (up and down) NJL model with SU{2) x SU (2) 
chiral symmetry and SU(N) color symmetry, with scalar and pseudoscalar couplings on 
the four-dimensional hypercubic lattice. We consider both naive and Wilson fermions and 
we study the model using a large N expansion as well as a Hybrid Monte Carlo (HMC) 
1^ numerical simulation. In this paper we will only present the main results. A detailed 
version of this work will appear elsewhere [||. 

The Lagrangian density in Minkowski space and in continuum notation is: 

C = ^(z^ - mo)^ + ^ [(^^)' + (^^75r^)'] . (1) 

In the above expression all indices have been suppressed. The fermionic field ^ is a flavor 
SU (2) doublet and a color SU{N) A^-column vector. The Lagrangian is diagonal in color, 
in contrast with the full QCD Lagrangian which is diagonal in flavor, r = {ti,T2,t^} are 
the three isospin Pauli matrices, ^ = j^df^, and mo is the bare quark mass (if mo 7^ 
the chiral symmetry is explicitly broken). To obtain a Lagrangian that is quadratic in 
the fermionic fields we introduce the scalar auxiliary field a and the three pseudoscalar 
auxiliary fields tt = {vti, 7r2, 713} [Q. 

Because of the chiral couplings the fermionic determinant has a phase. This phase is 
related to the chiral anomaly and the Wess-Zumino term and we will not consider it in 
this work. Going to Euclidean space and appropriately discretizing the above Lagrangian 
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we obtain the model on the hypercubic lattice. On the lattice, as it is well known, we 
have species doubling. The doubling in the NJL model will be interpreted as a doubling of 
the color degrees of freedom. To treat this problem we add to the Lagrangian density an 
irrelevant operator (Wilson term) of the form y^^^^', a being the lattice spacing and r a 
constant. Wc consider the r = case where no effort is made to remove the doublers (naive 
fermions) and also the r 7^ case where the doubler masses are raised to the cutoff and 
the chiral symmetry is explicitly broken (Wilson fermions). With these considerations and 
after appropriate scaling of the fields and couplings, so that only dimensionless quantities 
appear, we obtain: 

Z = J [d^d^dad7r]e-^ 

(N/2 ) 

x,y [i=i ^ ^ ) 

^xy = \Y.[i^i^~^)^x+n,y-{'yn + r)Sx-^,,y] + {4r + mQ + ax + i-f5'^x-T)Sxy (2) 

with 7^ hermitian, Uf — 2 the number of fiavors and /3i = 2n^- r — and mo — the 
symmetry is not explicitly broken. In that case at /3i = the model undergoes a second 
order phase transition from a symmetric phase with (a) — 0, massive tt and a fields and 
massless quarks, to a phase with spontaneously broken symmetry with (cr) 7^ 0, masslcss 
pion fields (Goldstone particles) and massive a and quark fields with dynamically generated 
quark mass equal to {a). In this work we will always stay in the broken phase. 

We study the above action using large N techniques and obtain analytical results both 
on finite and infinite volumes. The infinite volume results arc obtained sufficiently close 
to the continuum limit using asymptotic expansions. We also study this action for = 2 
using an HMC numerical simulation with Conjugate Gradient and leap-frog algorithms. 
The seven important results that stem from our analysis are presented below. 

I 

For naive fermions, at large N and with pion mass = 140 MeV, we calculate the a 
mass (Mcr), the a width (r^) and the quark mass (Mg) in physical units as functions of the 
a mass rrifj in lattice units. The results arc given in Figure 1 in units of the pion decay con- 
stant = 93 MeV. By setting Mq = 310 MeV we find = 726 MeV, F^ = 135 MeV, 
and A = vr/a = 1150 MeV (dotted vertical line). is consistent with phenomcnological 
expectations and A is consistent with the expectation that the cutoff should be close and 
below the mass of the lightest glueball (1550 MeV). The width however is underestimated. 
The reason is traced to the fact that to leading order in large A^ the width receives contri- 
butions only from the quark bubble and not from the pion bubble because the pion bubble 
is of order 1/A^. Because the phase space available for the a to decay to two quarks is much 
smaller than the phase space to decay to two pions, the pion loop contribution, although 
of order 1/A^, is probably more important than the quark loop contribution. 

II 

The above result is relevant not only for the low energy QCD but also for the Higgs 
sector. It is well known that there is an equivalence between the a-ir sector of QCD with 
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the scalar sector of the Minimal Standard Model. In the former the scale is set by the pion 
decay constant {F^, = 93 MeV) and in the latter by the weak scale (F^ = 246 GeV). As 
mentioned in I we find that in accordance with phenomenological expectations in the ct-tt 
QCD sector M^/F^ ~ 8 but in the Higgs sector all previous analysis predicts a triviality 
bound of the Higgs mass with Mo-/-F^ ^ 2.8 (see for example [0]). In the past this 
has been a reason for concern since it could imply that the Higgs mass bound may be 
underestimated. Our analysis suggests that this apparent discrepancy appears because the 
Higgs mass bounds are traditionally obtained for m„ < 0.5 while the M^/F^^ ^ 8 ratio is 
obtained for ~ 2 and it should therefore be accompanied by large deviations from the 
low energy behavior. Nevertheless, this is only a suggestion since we have not calculated 
deviations from the low energy behavior of a physical process that would enable us to make 
exact statements. However, the value of the width serves as an indication of the size of such 
deviations. In a way, departure from low energy behavior will be signaled by an increasing 
width of the a to two quark decay. At rria- ~ 2 the width is already fairly large. 




Figure 1. r = 0, -jr- = 93 MeV ^^'^ number of colors is set equal to three. From top 
to bottom the lines correspond to the a mass, quark mass and a width calculated at large 
and to leading order in mo-. The vertical line denotes the point where ^ = 93° j^e • 

III 

If the Higgs sector is the low energy effective field theory of a NJL model with = 3, 
Uf = 2 and = (the pions should now be viewed as the would-be Goldstone Bosons), 
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we obtain for the larger values of a figure that is almost exactly the same as Figure 1 
with (Mg.), {T„) and (Mg) measured in units of F^^ = 246 GeV. Strictly speaking such a 
figure does not contain any physical predictions but there is a very interesting observation 
that can be made. 

If we set the fermion mass to ^ ~ as is the case for the low energy sector of QCD, 
then the Higgs mass will be Mo- = 1915 GcV. This point will correspond to nifj = 2 as 
denoted by the dotted vertical line in Figure 1. As discussed in II, at rricr = 2 one would 
expect very large deviations from the low energy behavior of scattering cross sections. 
This suggests that as the CM energy is turned up first the deviations from the low energy 
behavior will become large, signaling the onset of new physics, and later the Higgs particle 
would be observed. 
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Figure 2. The k = 1/(8 + 2mo) plane for r = 1, = 2. The solid lines are constant 
lines. From top to bottom they correspond to m^r = 0,0.1,0.2,0.3. The dotted lines 
are constant niq lines. From right to left they correspond to niq — —0.2, —0.1, 0.0, 0.1, 0.2, 
0.3. The ruq — 0, ttItt = point is located at /^i^^^^^j = 0.3416, HcMrai — 0.3994. 

IV 

With Wilson fermions we obtain at large N analytical expressions of the pion mass (rriTr) 
and quark mass (mg) in lattice units as functions of the bare parameters of the model. We 
are then able to make exact statements regarding the approach to the continuum and chiral 
limits. For small nig and m^r the constant niq and lines including the nig — and m^^ — 
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lines are presented in Figure 2. The point Kchirai, where the ruq = and m,r = 

hnes intersect, is denoted by a circle and corresponds to the true chiral limit. 

Prom this figure we see that if for a fixed (3i we were to change k — ^^.^^^ from smaller to 
larger values (as is often done in QCD with dynamical Wilson fermions) then if Pi < 
we would reach the = limit before we reach the continuum limit. On the other hand 
if /5i > /5i , - , we would reach the continuum limit before we reach the m^r = limit. As 
mentioned above there is only one point in the k, plane where we can obtain a true 
chiral limit. This may provide an insight on how the retrieval of the chiral limit is achieved 
in QCD. 




1.5 2.0 2.5 3.0 

Figure 3. r = 0, itlq = 0, N — 2. The diamonds are the values of rriq and the crosses 
are the values of m'^ from the numerical simulation for Ix = 8, It = 16. The solid lines 
are the large N numbers with the zero mode included. From right to left they correspond 
to (L, = 8,Lt = 16), (L, = 16, Li = 16), (L, = 32, = 32), (L, = 64, = 64). The 
dotted hnes are the large N numbers with the zero mode excluded. Prom right to left they 
correspond to {L^ — 16, Lt — 16), {L^ — S,Lt — 16). The solid vertical line denotes the 
infinite volume Pi^ from large N. For Pi < Pi^ the model is in the broken phase. 

In Figure 2 the constant m^r lines were drawn under the assumption that <C Arrig 
which can not be satisfied on the niq — line, except on the one point where it intersects 
the rrij, = line. The constant m,r lines are therefore valid only in the regions where 
ml -C 4m^. 
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V 

At large N and for Wilson fermions the a particle has mass proportional to the cutoff. 
Our analysis traces this fact to two related reasons. First, although the Wilson term 
has raised the masses of the doublers to the cutoff, it has not decoupled them from the 
theory. Through vacuum polarization these contribute to the a self energy and raise its 
mass. Second, although the Wilson term has not altered the low frequency behavior of 
the propagating quark, it has however altered its high frequency behavior. Again through 
vacuum polarization the high frequency modes contribute to the a self energy and also 
raise its mass. Such a phenomenon may also be responsible for the difficulty in observing 
the (7 particle in numerical simulations of QCD with dynamical Wilson fermions. 
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Figure 4. a) The a propagator in momentum space for 10 small momenta with r = 0, 

rriQ = 0, (o") = 0.4840, N = 2, = 16, Lt = 16. The crosses are the values from the 
numerical simulation. The diamonds are the large N results, b) Same as in a but for 
(cr) = 0.35. c) Same as in a but for the pion propagator, d) Same as in b but for the pion 
propagator. 



VI 

The numerical and the leading order large N results are in good agreement, indicating 
that the 1/A^ corrections are small for the quantities we were able to measure. The first 
indication that the model at A'" = 2 has same type of behavior as at A^ = oo appears in 
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Figure 3. There we plot rrig and = {—N < > /(2/5i)} + mo vs. Pi as determined 
from the numerical simulation. At large N the gap equation gives rUg = m'^. This relation 
is satisfied nicely. Furthermore, the agreement with large of dynamically determined 
quantities vs. bare quantities is also quite good and helps us to get oriented in the bare 
parameter space. Figure 3 is an example of this. Some more results that demonstrate this 
agreement are presented in 
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Figure 5. r = 1, = 2, L3; = 8, = 16. The crosses are the MC data. The solid line is 
the large A^ prediction on the same size lattice. 

The important comparison with large A^ that will help us get a feel for the size of the 
1/A^ corrections comes from comparisons of dynamically determined quantities vs. other 
dynamically determined quantities. In particular we exchange one of the bare parameters 
for (a). In Figure 4 we present the a and tt propagators in momentum space for ten small 
momenta and r = 0. It is from these figures that we would have to extract the pion wave 
function renormalization constant Z^^. As it can be seen the large A^ prediction for the 
same (cr) as the one measured in the simulation is in good agreement with the numerical 
results. The large A^ predictions in Figures 4a, 4b, 4c, and 4d "fit" the numerical results 
with per degree of freedom 0.38, 0.67, 0.32, 0.42 respectively. This means that the pion 
decay constant in lattice units f-„ = (cr)y^ as a function of {a) has very small 1/A^ 
corrections. Another dynamically determined quantity that agrees very well with the large 
A^ prediction when plotted vs. (cr) is the pion mass m.^-. This plot is shown in Figure 5. 
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This figure suggests that the corrections to m.^ are very small. 

At large N we find that on a finite lattice can become zero only if (3i becomes infinite 
(free field). On the other hand we find that on a finite lattice can be made to vanish by 
adjusting the bare parameter tuq. Indeed we can see from Figure 5 that can be made 
very small. In fact it turns out that the value of mo where this happens is predicted by 
large quite well. 

As it has already been discussed in V, for r 7^ the a mass is of order cutoff and therefore 
very heavy to be able to measure from the decay of the a — a correlation function. However 
for r = one would expect to be able to measure mo-. As we will discuss in VII below 
this is not possible with the lattice sizes accessible to us. This is unfortunate since mo- is 
another very important quantity. However, we expect the size of the 1/A^ corrections of 
mo- to be similar to the ones of m^^ and therefore very small. Also, we did not measure the 
cr width in the numerical simulation, but, as discussed in I, we expect the 1/A^ corrections 
to be large. 

VII 

On finite lattices the zero momentum modes of the quarks affect the inversion speed of 
the CG algorithm and introduce finite size effects. To leading order at large A^ the matrix 
M^M of eq. ^ is diagonal in momentum, spin, flavor, and color spaces. The smallest 
eigenvalue of this matrix is m^ and corresponds to the p = matrix element. This in turn 
corresponds to the zero momentum mode of the quarks which from now on we will simply 
call "zero modes". 

On a finite lattice the role of the zero modes is crucial in the inversion of M^M. For 
small rriq, the condition number of the matrix is 4/m^ for r = and 64r/m^ for r = 1 and 
it is clear that it depends strongly on the presence of the zero modes. A large condition 
number will make the inversion of M'^M very slow. Furthermore the spacing of the smallest 
eigenvalues behaves like and for larger lattices the inversion times will rapidly get 

worse. An important observation can be made by noticing the dependence of the condition 
number on r. This suggests that performing the simulation with smaller r will yield a quite 
faster inversion. It is possible that this may also be the case for QCD. 

But the unwelcomed effect of the zero modes on a finite lattice is not limited to large 
inversion times. Because on a finite lattice their effects are not suppressed by the measure 
but instead by an inverse volume factor, it turns out that in certain cases they severely 
obscure the physics. By simply looking at the numerical results in Figure 3 we would 
not only be unable to estimate the critical point but also we would be unable to see any 
indication of a phase transition. The large A^ result on the same size lattice also has the 
same problems. As we increase the lattice size in the large A^ calculation (solid lines from 
top to bottom) we see that a picture of an order parameter slowly materializes. At = 64, 

= 64 a fairly good prediction of the large A^ infinite volume critical point is achieved. 
If we now do the same large A^ calculation but neglect from the momentum sums the zero 
modes, we obtain as a result the two dotted lines for = 8, Lt = 16 and = 16, Lt = 16 
(from left to right). We see that neglecting the zero modes on an = 8, Lt = 16 lattice 
gives very similar results as the ones obtained on a L^. = 64, Lt = 64 lattice with the zero 
modes included. 

As mentioned earlier although one would expect to be able to measure the a mass in the 
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r = case we were not able to do so. Large N provides an explanation of this unexpected 
problem. If we plot the real part of the inverse a propagator for a finite lattice and set the 
external momentum to g = {im„, 0, 0, 0} we can obtain the sigma mass at the zero of this 
function. We find that because of the presence of a discontinuity we do not obtain a root 
at all. A root will eventually be obtained but not until rricr becomes very heavy [Q. The 
presence of this discontinuity is again due to the zero modes. 

There are some very interesting issues relevant to lattice work that have not been 
considered in this paper. It would be important to calculate the three and four point vertices 
and therefore be able to calculate scattering amplitudes and their departure from low energy 
behavior as well as the corrections to the width. It would also be interesting to study 
the NJL model at finite temperature and investigate the finite temperature transition in 
connection with the approach to the chiral limit. Finally it would be important to include 
vector meson couplings (see, for example, 0,0) and confirm that the vector meson masses 
scale appropriately and do not become of the order cutoff as the a particle does. 
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